function [Fval,Model] = firstOPoly_Poiss(theta,Xvec,Data)
%FirstOPoly_Norm Sum squared error for 1st order polynomial

a=theta(1);
b=theta(2);

Model=a+b*Xvec;
Fval=sum(Model-Data.*log(Model));

end

